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1 Introduction 


In the light of the recent progress in the four-loop calculations [1, 2] one of the topical 
tasks is the construction of e-expansion of one-, two-, and three-loop master-integrals of the 
existing packages like ONSHELL2 [3] or three-loop packages [4, 5, 6, 7] up to level sufficient 
for getting finite four-loop corrections in physical quantities. The all-order e-expansion of 
the one-loop propagator diagram with arbitrary masses and external momentum, and the 
two-loop bubble diagram with arbitrary masses was constructed in [8, 9]. The results are 
expressible in terms of Nielsen polylogarithms [10] only. The construction of the e-expansion 
for diagrams relating to QED/QCD problems was investigated in [11, 12], The finite parts of 
the three-loop bubble integrals were presented in [13, 14, 15]. The higher order e-expansion 
for some of the master-integrals from these packages was calculated in our previous papers 1 
[9, 16]. Independently of analytical calculations, the numerical approach to the evaluation of 
the single scale diagrams has been developed in recent years [17]. Based on this technique, 
the high-precision numerical values of higher order coefficients of the ^-expansion for four- 
loop bubble integrals [18] and three-loop propagator type integrals on mass shell [19] have 
been calculated recently. 

The present paper is devoted to the analytical calculation of the higher order terms of 
the e-expansion for the scalar integrals shown in Fig. 1. At M = m these integrals enter 
in packages ONSHELL2 [3] (Vim) and three-loop packages [4, 5] (D 5 ). The numerical 
values of higher order coefficients of the e-expansion for diagram Vmi and D 5 were cal¬ 
culated by Laporta in [17]. Another problem under consideration is further investigation 
and development of the “sixth root of unity” approach proposed by Broadhurst in [13] and 
developed in [14, 16, 9, 20, 18]. The main idea is that transcendental numbers occurring in 
the e-expansion of single scale diagrams (diagrams with only one mass scale) are defined by 
the values of massive cuts. For the diagrams with zero-, one-, two-, and three- massive 
cuts, the set of transcendental numbers are related to generalized log-sine functions [21] or 
their generalization [9, 20] of special values of arguments {7r, 7 t/2, 7t/3, 27t/ 3}. For diagrams 
with four massive cuts new constants, associated with the elliptic function [18] appeared. 

We work in the dimensional regularization [22] with space-time dimension n — 4 — 2e. In 
our normalization each loop is divided by (47r) 2_e r(l + e). We also use the following short 
notation for the auxiliary integrals appearing in our calculation: 


J, 


771177127713 


V r 


771177127713 


d n kid n k 2 


B 0 (mi,m 2 ,m 3 ) = j 


[(h-p) 2 +ml] [kl +m|] [(/q - k 2 ) 2 +ml] 
d n k\d n k 2 

[k\+ml][k%+ml][(k 1 -k 2 ) 2 +ml] 
d n ki 


[(ki—p) 2 +mi][k 2 +m 2 \ 


p 2 =—m | 


Where is a typo in Eq. (4.10) of [9]: the coefficient before 7 tLs 4 (|) should be instead of We are 
grateful to Y. Schroder and A. Vuorinen for correspondence. 
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2 Vr 


mmmM 


This integral enters in two-loop relation between pole and MS masses of heavy particles like 
t-quark or Higgs boson within SM [23, 24], The integrals with an arbitrary set of indices, 


b mmmM ( O, , (T, A ) / / / ; \ 

J ((p — ki) 


d n kid n k 2 

2 + m 2 Y((ki — k ; 2 ) 2 + m 2 ) a (k 2 + m 2 )P(k 2 + M 2 ) x 


where the external moment belongs to the mass shell, p 2 = —m 2 , can be reduced to one 
master-integral with indices (1,1,1,1). For the calculation of the diagram, we use the differ¬ 
ential equation technique [25] 


i YnmmM (1 j 1 j 1 j 1) b mmmM (1 j 1 j 1 j 1) 


- V Mm m + 2 A 0 (m)B 0 (M, m, m 


1 / n 


1 / 3n 


m 2 V 8 


’ ’ \ 4m 2 V 2 

—---'j 

M 2 M 2 — Am 2 ) 


1 (n - 3) 
M 2 + M 2 - 4m 2 


1 1 
M 2 _ M 2 - 4m 2 


The analytical solution of this equation up to part linear in £ was presented in [24]. Equation 
(2.3) has a very simple form in the framework of the geometrical approach [26]. Let us 
introduce a new variable 


cos 6 = -, M < 2 m 

2m 


2 



Taking into account all order e-expansion of the one-loop propagator and two-loop bubble 
integrals [27, 8] we rewrite the original differential equation (2.3) 


1 d 
sin 6 cos 6 d6 


V r 


mmmM 


- 2 


£ 1 -2e 

cos 2 9 sin 2 9 
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V r 


mmmM 


2 cos 2 0 sin 2 0(1 — 2e) [ e 2 


3 — 6 cos 2 9 + 2(4 cos 2 d) 


l—£ 


20 00 (2 fV 

-|-(4 cos 9 sin d) 4 ~ 2£ + (4cosdsiiid) 4 ~ 2c E [4Lsj +2 (2 9) — 3Lsj +2 (4 9)\ 


=o (j + 1)! 


j =o 


J 9 _ v 

° mmm ^ uo 


= o. 


m 2 2 sin 2 9 cos 2 9 
For a modified function, V (9) defined as 

Vm.mmM = (sind) 2 ~ 4e (cosd)~ 2£ (m 2 )~ 2£ V (9) , 
the solution has the following form: 


(2.5) 


( 2 . 6 ) 


2(1 - 2e) 2 C(d) = -(2 - 3e)(l - 2e)- 


J >7 


(m 


2U-2e 


sind) 4£ (cosd) 
2 sin 2 9 


(1 - 3e)Isc 4i2 (9) 


2el 


3(1 - <r)Isc 4 , 2 rn + Ihnhiw: _ 


.,_ 2 “ ( 2 e ) J -‘ 

+ 4 HHi— 

l=o 


cos 0 

1^21 [4Ls J+ i (2d) - 3Lsj + i (4d)] 


"2~ 2e E 


k=0 


M 

fc! 


[ (sin d) 
,ln fc+J ' +1 (2 sind) 


+ 3eLsfc+i (2d) Lsj+i (4d) 


_22-2e £ ^ (2 £ ) A f d Toi„fc 


k =0 


k\ 


k + j + 1 

dtp 12 ln ft (2 sin 0) Lsj + i (20) + 3Lsfc +4 (20) [In (2 sin 0) + In (2 cos 0)] J 


+3 • 4 1 £ V ^~rr~ / dtp——— ln^ (2 sin 0) [In (2 sin 0) T In (2 cos 0)P \ , 

k=o kl J sin ^ J 


(2.7) 


where Lsj (d) is a log-sine function [21] and we have introduced a new function Isc a ,b (d) 

r 9 (sin 0) ae (cos <p) be 


Isc a , b (d) = J dtp- 


sin 0 cos 0 

where a, d are integer numbers. The e-expansion of this function is 

i - (E* 

Isc a , b (d) = In (tand) + - E 7 i~—pr. (a 7 hT +1 (sin 2 9) - b l ln t+1 (cos 2 9)) 

1 i=i \ z + 

1 OO / 

~T) E ( _ b ) (a*Si,i(cos 2 d) - d*S M (sin 2 d)) 


( 2 . 8 ) 
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£ M 


_ £ \ (i+j) r ‘ (—I)’’ Ill' (sin 2 




(sin 2 e) 


-E 


(—1 ) r ln r (cos 2 9) 


Sj+l— r,i (COS 9^j , 


where S a ,b(9) is the generalized (Nielsen) polylogarithms [10]. The last integral in Eq. (2.7) 
can be explicitly integrated 

/ 0 f*OS f) 1 

d6 —— ln a (sin 6) =-ln a+1 sin 9 , 

sm 9 v ’ 1 + a 


r e cos9 a b 6 !(-l) b “ a!(-l)* 

dQ --In (smd In (cosd -—— > -- 

sin 9 v ’ v ; 2 1 + a + 6 ^ (a-i)\ 


6 ) = E §^77 ( 2 ln ( sin 0 ))° 1 S *+^ ( sin2 0 ) ( 2 - 10 ) 


The solution for the original integral can be represented as 


V mmm M = (sill 44 (cos 0) 2s (m 2 ) 2 ‘\V(0)+ £ £ 


( 2 . 11 ) 


where the coefficients Vj should be defined from the boundary condition. We choose a 
particular value of the diagram at M = 0 which corresponds to 9 — n/2. In this case the 
diagram reduces to J mmm and the product of one-loop vacuum diagrams 


KtrnimM | M =0 


3n - 8 Jmmm _ 3(n - 2) 2 A 0 (m 2 ) 

4(4 — n) m? 8{n — 4)(n — 3) m 2 


Using the results of [11, 28] we get the following values for the first several coefficients: 

v-! = 12 In 2 , v 0 = 48 ln 2 + 6(2 - 36 In 2 2 , 

Vl = 144 ln 2 + 24C 2 - 5Cs - 12 ( 2 In 2 - 144 In 2 2 + 72 In 3 2 , 
v 2 = 3841n2 + 72C 2 -20C 3 -29C 4 -48C 2 ln2 + 24C 2 ln 2 2 
-432In 2 2 + 288In 3 2-106In 4 2+48Li 4 (±) , 


v 3 = 9601n2 + 192C 2 -60C3-H6C 4 -^C5-144C 2 ln2 + 96C 2 ln 2 2 

-48C 2 In 3 2 + 306C 4 In 2 -1152 In 2 2 + 864 In 3 2 -424 In 4 2 + ^ In 5 2 


+192Li 4 (|) + 288Li 5 (| 


( 2 . 12 ) 


v 4 can also be calculated by using the results of [28]. 

From Eq. (2.7) it is evident, that the coefficients of the e-expansion of the diagram U m2 
can be parametrized by introducing a new class of functions LsLsc^j ( 9 ) defined as 


LsLsc^jj ( 9 ) = / d(j) Ls fc+ i ( 0 ) In* 1 2 sin — In - 7 1 2 cos ■ 


(2.13) 


where k,i,j are integer numbers, k> 0 and i, j > 1. Some properties of these functions are 
collected in Appendix B. 
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Let us write the ^-expansion of the diagram in following form: 




2\-2e 


+ -S(6)+F(0) +eE(9) +£ 2 N(9) +£ 3 N n (9 ) +£ A N nn (9) + 0(£ 5 
l£ z £ 


(2.14) 


The functions S(9),F(9 ) and E(9) were calculated in [24] (see Eq.(3.27)). The analytical 
results for N(9) for arbitrary values of the angle 9 read 


N(9) = l F(9)- E(9) _ 4Q 


sin 2 9 


sin 9 


4Ls 3 (2 9 ) - 3Ls 3 (40) 


|^E(-l) , l^4- j (2^«3Ls 4 _U40) 


3=0 


lie 


1 


sin 2 9 


f + f cos 2 9L A + 24C 2 hi 2 - 14Cs - 12C 2 


-24$ (20) + 12Lsl 1} (20) - 3Lsf } (40) + 20 


1 2 


Ls 2 (20) 


12Ls 2 (20) Ls 2 (40) 


-160L e Ls 2 (20) - 320/ e Ls 2 (20) + 24C1 3 (vr - 20) [L e + 2 l e \ - 12C 2 Li 2 (sin 2 0) 

-80 2 L 2 - 320 2 / 2 - 32 9 2 L e l e + 36Cs^ - 24C 3 L e + 12C 2 L 2 + 48 ln2C 2 L 0 

+42C 3 In 2 - 31C 4 + 2 In 4 2 - 36C 2 In 2 2 + 48Li 4 (|) , (2.15) 


where 


Le — In (2 cos 0) , Iq — In (2 sin 0) , l m e — In (2 sinm0) , 


and 


$(0) = J d(j) Ls 2 (0) In [2 cos ■ 


LsLsci i 2 (0) , 


(2.16) 


is the function defined in Eq. (2.41) of [29]. The results for N n {9 ) and N nn {9 ) are sufficiently 
lengthy to be published here. It should be mentioned that in the order of e 3 the following 
combination appears: 


2LsLsci, 2j2 (20) + LsLsc o , 3 , 2 (20) + LsLsc o , 3 , 2 (7t — 20) + LsLscli. 3 (20) — LsLscip^ (it —29) .(2.17) 

To write results valid in other regions of the variable (M > 2m), a proper analytical con¬ 
tinuation of all expressions should be constructed. For generalized log-sine functions it is 
described in detail in [9, 31]. In terms of the variable 

y = e ia2Q , hi(—y — iaO) — \ny — ian, (2-18) 

the analytical continuation of all generalized log-sine can be expressed in terms of Nielsen 
polylogarithms, whereas for the function $(0) the result is expressible in terms of harmonic 
polylogarithms (see Eq.(3.4) in [29]) introduced by Remiddi and Vermaseren [32], In terms of 
conformal variable (2.18), the analytical continuation of the LsLsc functions can be written 
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in terms of harmonic polylogarithms (see the discussion in [9, 29, 33]). In particular, an 
analytical continuation of the LsLsc 1^3 (^-function produces integral / -ff-i,o,o,i(z)ypr- The 
result is relatively lengthy to be published here. The analytically continued results of the 
order e 1 checked by heavy mass expansion [34] with the help of the packages described in [35]. 

For a particular case M — m the integral V mrnrn M converts into the master-integral Vmi 
from the package ONSHELL2 [3, 30]. Its ^-expansion is (we present only new coefficients 
of the expansion) 


N (f) = Vim[e 2 ] = f - - 32^J + 12 In3 - 2 ^ In 2 3 + In’ 3 
-346 + I 6 C 2 In 3 - 46 hi 2 3 - 326 + fCs In 3 - 186^ + §673 In 3 + 96^ 

• Ls >(il + In 3 + f irL S2 (f) - 42 ^^ 


-84- 


V3 


LS; 


.1 Wal ln 2 3 _ |7rLs 2 (f) In 3 + f ^ ^ hi3 - 7 ttLs 3 (f) 


- 7- 


Ls. 


m 


V3 


+ f Lsf } (t) + T Ls 2 (f) 


and 


+246 In 2 + 18C 2 hi 2 2 - 186 In 2 In 3 + 9C 2 Li 2 (±) , 

N„ (f) = V ml [£ 3 ] = if - 80+ 32In3 - 6 ^ In 2 3 + |^ ln : > 3 - ^ !„■ 3 


(2.19) 


-54 < 2 ^ + 18<2^ In 3 - !< 2 ^ In 2 3 - 148(2 - 152< 3 - 157C 4 + ^( 2(3 + 
+48C 2 hi 3 - 16C 2 hi 2 3 + §6 hr 3 3 + 366^ - 96^ In 3 + 786 In 3 - f 6 In 2 3 

+636 In 3 — ^ 6 ^ ~ 224 Ls ^ + 84 Ls ! 3) In 3 - 14 LSa( ~^ ’ 9 " ’ 7 ' LS2f “ 


V+ ^ +3 

+ 167 rLs 2 - y 7 rLs 2 (j) In 3 + § 7 rLs 2 (j) In 2 3 + 67 r 6 Ls 2 ( 

LS3 


3Z ln 2 3+ l^nii ln 3 3 


6 v+ 


Ls E 


-126- 


M , i: L3 4t) 




V3 


In 3 — 


2 i 


<z) 


V+ 


In 2 3 - 14vrLs 3 (f) [ 2 -In 3] 


Ls 


+i 7rLs 4 (f) - f ■^ Ls 4 (f) _ 7 —t 4 - ln3 ] + 27Ls i 1} (t) [2 “ ln3 ] 


Ls. 




_i_ 28 

^ 3 


Ls 2 (f)] 2 [2 - In 3] - 4Ls 2 (§) Ls 3 (f) - f vrC 2 Ls 2 (f 


2 Vs 

+9Lsg ^ — f)t 5 T 36LsLsciy j3 + 1686 In 2 

-4 In 4 2 + 246 In 2 2 - 96Li 4 (±) - f In 5 2 - 63Cs In 2 2 

366 + fCs - I 8 C 2 In 2 - 186 In 3 


if- /i5 


Li 4 Q) ln 2 + Li 5 Q 
96Li 3 (|)-186S 1;2 (i 


+72Li 4 (|) ln 3 + 3 ln 3 In 4 2 - 726 In 2 ln 3 + 186 In 2 In 2 3 . 


( 2 . 20 ) 


Numerical values of the coefficients, Vmi[£ 2 ] and Vnnfe 3 ] are in full agreement with the 
proper results of [17]. 
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3 D 


mmmmM 


The differential equation for this integral is 


—M z 


d 


dM 2 


D 


mmmmM 


= 5 + 4(3 - n) 


m 


+ b n ( ~n - 4 


M 2 - Am 2 


4“ l / M mm 74o(nr 


M 2 - 4m 2 2 

4 — 2n 


- -n £ 


mmmmM 


M 2 - 4m 2 ’ 


(3.21) 


where Bn = .B/v(0, 0,1,1,1,1) was defined by Broadhurst in [11, 12]. Using again the angle 
variable defined in (2.4), we rewrite equation (3.21) in the following form: 


cos 6 d 
2 sin 6 dd 


D 


mmmmM 


-1 + 3e + 


1 — 2e 
sin 2 9 


D 


mmmmM 


B n (2 — 3e) 1 (m 2 ) 


2 \ 1 —3e 


m 


2 4 sin 2 6 sin 2 6 e{l — e)(l — 2e) \ e 2 


1 — 2 cos 2 9 + (4 cos 2 9 ) 


l—£ 


(2e) 


+ (4 cos 0 sill 0) -J2 jjUfjl l 2Ls 3+2 ( 2ff ) ~ Ls j +2 (49)11 , 

For auxiliary function D(9) defined as 

DmmmmM = (silld) 2 ^ (cOsd)~ 2£ (?7T 2 ) D{9) 

the solution is 

(l-e)(l-2. efD{9) = 

(1—e)(2 —3e)(l —2e) B N 


(3.22) 


(m 


2\2—3 e 


(1 - 3£)Isc 4i2 ( 6 ) 


sind) 4e (cosd) 2c 


2 sin 2 9 


£ 

-2 


(1 _ e)1sc 4 '2 (0) + (sm9)4E U, Sg)2g - 2 1 ~ 2s (slll9)4 ‘ 


2 sin 2 9 


sin 2 9 


3-4sl v_(2 e) j 

2-^ (A 


£P)(j + iy l (sin 9) 


COS f) 

-1 *5r |2 Ls, +2 (20) - Ls_,+2 (46*)] 


,-2 e f. ( 2 £) 

h H 


" 4 ln ^ J+2 ( 2sing ) + eLsfc+1 ( 2 9 ) Lsj +2 (40) 


k + j + 2 


-2 2 " 2 V]T 


(2f) 

/.-! 




d(j) Lsfc + i (20) [In (2 sin 0) + In (2 cos 0)] J+1 + ln fc (2 sin 0) Lsj +2 (20) 


+2 2_2e £ 


k =0 

00 (2e) fc r° 


k =0 


A:! 


d(j) 


1 ° ln fc (2 sin 0) [In (2 sin 0) + In (2 cos 0)] 2+1 1 , 


sm 1 


(3.23) 


where the last integral can be explicitly integrated with the help of Eq. (2.10) and Isc a) 6 (9) 
is defined in (2.9). The result for the diagram is 


DmmmmM = (sin 9) £ (cos 9) £ (m 2 ) 1 3£ j D(9) + £Jd j | > 


(3.24) 
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where the coefficients dj are defined from the value of the original diagrams at M = 0, which 
corresponds to 6 = 7t/2: 


m 2 D 


mmmmM 


M =0 




n — 2 
n — 4 


— Bn 


3 n — 8 
4 (n — 4) 


(3.25) 


The explicit values of the first several coefficients are 


d_ 3 = §, d_ 2 = y — 8 In 2 , d_ 4 = f - f & - 40 ln2W 8 In 2 2 % 

d 0 = + f C 3 -- f C 2 - 136 In 2 - 40 In 2 2 + If In 3 2 . 

d x = - i n 4 2 + 8 §o ln 3 2 + I 6 C 2 In 2 2 - 136 In 2 2 - 392 In 2 

+ ffs + l -fU ~ 64Li 4 (|) + 172 , 
d 2 = 428 + ln 5 2 _ 2000 ln 4 2 + 2^2 ln 3 2 _ 392 ln 2 2 _ 1032 ln 2 

-32C, In 3 2 + 8 OC 2 hi 2 2 - 408C 4 In 2 - ^C 2 + Cs + 4 

+f C 2 C 3 + 422^5 - 320Li 4 (1) - 384Li 5 (§) , (3.26) 

The next coefficient d 3 contains the U 51 constant [12] which could be rewritten in terms of 
the generalized log-sine function of the argument 7 t /2 [9]. Collecting all the results we get 
the first several coefficients of the e-expansion of the diagram 


(m 2 ) _1+3e (l — e)(l — 2e) 2 D. 


mmmmM \ 


( 0 ) = - J3 {1 + 2 cos2 9 } ~ ^2 I 1 " 12 cos2 6Le } 


3e 2 


-- {l + 4cos 2 0Ljj - 2sin(20) [2Ls 2 (20) - Ls 2 (40)]} 


16 

-6 H-cos 2 OLl — 4 sin 2 0 

3 


4C1 3 (20) - Cl 3 (40) - -C 3 


+4sin(20) [2Ls 3 (20) - Ls 3 (40)} - 8 sin (2 0)l 2e [2Ls 2 (20) - Ls 2 (40)} 

{ 56 

— 18 + -^-(3 + 8 sin 2 0 [2lg + Lg} [4C1 3 (20) — Cl 3 (40)} 

+4 sin 2 0 [8F(20) - 4Lsf ) (20) + Lsf 5 (40)" 

+16 sin 2 0 0 [2Ls 3 (20)-Ls 3 (40)]-Ls 2 (20) [2Ls 2 (20)-Ls 2 (40)] + ( 3 (2L e -3lg) 


+4 sin 2 0 


17G - 16Li 4 f 1) + 4 ( 2 In 2 2 - 14( 3 ln 2 - - In 4 2 
w 3 


■ cos 2 OLl 


+ 8 sin( 2 0)4 [2Ls 2 (20) - Ls 2 (40)] - 8sin(20)/ 2e [2Ls 3 (20) - Ls 3 (40)] 
+^ sin(20) [2Ls 4 (20) - Ls 4 (40)]) + 0(e 2 ) , 


(3.27) 


where l m g and Lg are defined in (2.16). The analytical results for the next two coefficients are 
expressible in terms of the LsLsc-functions. In the order e 2 the previous combination (2.17) 
of LsLsc functions is included. The divergent parts of this integral were calculated previously 
in [36]. In the regions of the variable (M > 2m) the proper analytical continuation of all 



expressions can be constructed. It complctly coincides with the previous case V mmm M • The 
result is relatively lengthy and therefore will not be presented here. 

For the case of equal masses M = m (the master-integral D 5 ) we get 


[m 


2\-l+3e 


(1 - e)(l - 2 e) 2 D 5 = 1 - 6 


3e 2 

LS; 


,Ls 2 (?) 


-6+ 6C 2 ^ +6(3 + 18^5^- 


12 


+3 

Ls 2 f? 




22 In 3 


e{ 18 + 12C 3 ^ + 6C 2 ^1h 3 - f C 3 + 26C 3 ln3 - 63C 4 - 6 


y/3 


56, 
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Ls 2 f? 




^2 hr 3 
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+18 y > In 3 - 12vrLs 3 (f) + 18Lsf } (f) + 8 [Ls 2 (f) 
+£ 2 |-54+3C 2 ^ln 2 3 + 12C 3 ^ln3 + 56C3-136C4-^C2C3-^C5 


2 LS 4 

- 12 
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Ls 2 
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(t) 


+26C 3 In 2 3+f C 4 ^ - 24C 4 In 3 - 2 3 + 9 ' hr 3 -12 


Ls, 


+3 


It) 


V3 


In 3 


+16 Ls 2 


! 2 Ls 5 

In 3 + 6 


ffi 


38 


7rC 2 Ls 2 (|)-12Ls| 1} (y) + §X5 


+3 


-24vrLs 3 (f) ln3-f ttLs 4 (f ) + 12vrLs 4 (f)+36Lsi 1} (f) ln3-48LsLsci, li3 (f) 
-32( 2 In 2 2 + f In 4 2 + 128Li 4 (±) 

—64(^ 2 In 3 2+84^3 In 2 2+y In 5 2+384 [ln2Li 4 (±) In2 + Li 5 (|)' 

+24^2 In 2 2 In 3 —42£ 3 Li 2 (|)-96Li 4 (|) In 3-4 In 4 2 In 31 . (3.28) 


As a non-trivial check of these results we established full agreement between the numerical 
values for the coefficients of e-expansion of the Lb-integral with the proper Laporta results 
[17], 


4 Conclusion 

In this paper, the higher order e-expansion of the diagrams shown in Fig. (1) have been 
constructed (2.7), (3.23). The coefficients of the e-expansion are parametrized in terms of 
generalized log-sin functions and LsLsc-functions described in Appendix B. At M = m these 
integrals enter in FORM [38] based packages for calculation of two-loop on-shell self-energy 
diagrams [3] (Vim) and three-loop vacuum integrals [4, 5] (L> 5 ). The numerical values of the 
calculated coefficients (2.19),(2.20) and (3.28) coincide with the results presented in [17]. 

We shown that the basis of transcendental numbers for single scale diagrams with two- 
and three-massive cuts contains new elements in addition to the odd/even basis of weight 5 
constructed in [14, 9]. Some of these elements are the product of the lowest weight elements 
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of the ” sixth root of unity” basis: 

C2 x {in 2 In 3, In 2 In 2 3, Li 3 (±) , Si, 2 (j) } , Li 2 (|) x {C2, Cs, C2 ln2, ( 2 In 3} , 

Li 4 (i) x {In 2, In 3} , In 3 In 4 2. (4.29) 

Only one extra term should be added to the weight 5 basis. This new constant can be related 
to a LsLsc^i^-function of the argument |. The results of analytical calculation of three-loop 
master-integrals, evaluated in [9] and in the present paper are available on 

http : //theor.j inr.ru/~kalmykov/three — loop/master.uu 


Acknowledgements. I am grateful to A. Davydychev, F. Jegerlehner, and O.V. Tarasov 
for useful discussions. I thank Y. Schroder and A. Vuorinen for correspondence concerning 
the typo in expression for the e-part of D 4 and for information on the results of their 
calculation before publication. I would like to thank S. Mikhailov for his interest in work, 
and G.Sandukovskaya for careful reading of the manuscript. This research was supported 
in part by RFBR grant # 04-02-17192 and the Heisenberg-Landau Programme and NATO 
Grant PST.CLG.980342. 

A Auxiliary integrals 

Below we present the values of some integrals appearing in Sec. 2. 
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B Auxiliary function 

Let us describe here some properties of the new function LsLsc^jj ( 9 ) defined as 

LsLsCfc i j (9) — [ d<p Lsfc+i (0) In 1-1 2sin — In- 7-1 2cos— , (B.l) 

Jo 2 2 

where k,i,j are integer numbers, k > 0 and i,j > 1, and 9 is an arbitrary real number. 
These functions have appeared in higher-order epsilon expansion of two- and three-loop 
Feynman diagrams investigated in this paper (see Eqs.(2.7) and (3.23)). LsLsc-functions 
are also related to the e-expansion of hypergeometric functions. As was shown in [20] the 
multiple inverse binomial sums ^Si(n— 1) and ^S\(2n— 1) , where S a (n) 

is the harmonic sum, are expressible in terms of LsLsco^j (9), LsLsci^i (9) and LsLscipj (9). 
Moreover, the elements X 5 and X 5 contributing to odd/even basis of a set of transcendental 
numbers [9] are related to LsLsci )3> i and LsLsci j3j i respectively (see Eq. (15) in 
[20]). LsLsci^i (9) also appears in higher order e-expansion of the one-loop triangle diagram 
(see Eq. (4.5) in [29]). For lower values of the parameters, function (B.l) reduces to the 
generalized log-sine functions. In particular, 

• for k = 0 and i — 1 or j — 1 

LsLsc 0 ,i,i ( 9 ) = Ls& ( 9 ) , (B.2) 

LsLscoqj ( 9 ) = Ls^ (i t — 9) — Ls^ (n) — tt [Lsj (n — 9) — Lsj (7r)] . (B.3) 

• for k — 0, i = j — 1 

LsLsc 0 , 2,2 (0) = |Lsi 1} (29) - iLsi 1} (0) ~ j [Lsf } (tt -9)- Ls^ (tt)' 

+ |7T [Ls 3 ( 7 T - 9) - Ls 3 (7t)] . (B.4) 

• for an arbitrary k and i = j = 1 

LsLsc M)1 (9) = 9Ls k+1 (9) - Lsj£ 2 (9) . (B.5) 

• for k — 1, i = 2 and j = 1 

LsLsc 1)2)1 (0) = -i [Ls 2 (9)} 2 . (B.6) 
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For k — 1, i — 1 and j — 2 we get a function, which was denoted as <&(#) in [29] (see 
Eq.(2.41) in [29]), 

LsLsc 1)1i2 {9) = $(0) . (B.7) 

The relation between the functions of the opposite arguments is 

LsLsc^j (—9) = LsLsc k4J (6) , (B.8) 

The symmetry relation for the new functions 

LsLsCfc^j (9) — LsLsc^jj (27 t — 0) + LsLsCfc i jj (27r) = 2Lsfc +i ( 7 r) [LsCjj(27r — 9) — Lsc ij -(27r)][B.9) 

includes the generalized log-sine-cosine integral introduced in [9] (see Appendix A.2). Its 
definition is 

Lsc i,j(9) = — J d<p ln i_1 2 sin ^ In- 7-1 2cos^ . (B.10) 

In particular, the following relation is valid: 

r6 

/ dcpLsCij(4>) = 9LsCij(9) — LsLsc 0 ,jj (9) . (B.ll) 

J 0 

At the same time, the relation (B.9) can be considered as definition of the generalized log- 
sine-cosine functions (Ls*, ( 77 ) is a normalization constant). For k = 0 there is an extra 
symmetry relation, 

LsLscoyj (9) — LsLscoy,i ( 7 r — 6 ) ) + LsLscoy,i (n) = — 7 r [Lsc^ (77 — 9) — Lsc^tt)] . (B.12) 

When i or j are equal to unity, the following relations are valid: 

LsLscfc^i ( 9 ) + LsLsCj-i^+i,! ( 9 ) = —Ls^.+i ( 9 ) Ls* ( 9 ) , (B.13) 

LsLsc mj (9) + LsLsCj.i^fc+i ( 7 r - 9) - LsLsc^i^+i ( 7 r) = Ls fc+ i (9) Lsj {it - 9) ,(B.14) 

The values of the function LsLsc^jj {9) of the argument 9 = 2tt can be extracted from 
relation (B.9). Using the symmetric properties of the LsCj j (#)-function, 

Lscij {—9) = -Lsc i tj {9) , 

LsCjj (27 t — 9) = LsCjj (27r) — Lsc,j (9) , 

Lscjj (27 t) = 2LsCij ( 77 ) , Lscj j ( 77 ) = Lsc^ ( 77 ) , (B.15) 

it is easy to get 

LsLsc k ,ij (27r) = -2Ls fc +i M LsCjj(7r) , (B.16) 
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where the values of Ls^+i (7r) and LsCjj {n) can be founded in Lewin’s book [21], Taking into 
account that Ls 2 (7r) = 0 we get for k = 1 


LsLsci^j (27t) = 0 . (B.17) 

Using relations (B.16), Eq. (B.9) can be rewritten as 

LsLsCfc^j (6) — LsLsCfc t jj (27r — 9) = — 2Lsfc + i (7r) [LsCjj(d) — Lsc;j(7t)] . (B.18) 

From relation (B.12) we get for 9 = n 

LsLsc 0 ,ij (7r)+LsLsc 0 j,i (7r) = 7rLsCjj(7r) . (B.19) 

For particular values of the LsLsc-functions of weight 5 at 9 — 7r, 7 t/ 3 and 6 = 2n/3 the 
PSLQ analysis [37] yields 


LsLsco, 3,2 (tt) — |C 2 hi 3 2 — In 5 2 —|C 3 hi 2 2 — 4 ln2Li 4 (|) + ^C .5 + |C 2 C 3 — 4Li 5 Q) 

= 1.1346108755961038391083966682727... , 

LsLsc 1)1;3 (tt) = § In 5 2 - 2( 2 hr 3 2 + f Cs hr 2 2 +12 In 2Li 4 (|) - R 2 C 3 - WCs + 12Li 5 (|) 
= 0.6301340120387385042438188494654... , 

LsLsc w (vr) = -j|ln 5 2 + §C 2 ln 3 2 - JC 3 In 2 2 - 41n2Li 4 (±) + - 4Li 5 (±) 

= 0.8874478318089418243337609880871... , 

LsLsc 4j2i2 ( 7 r) = In 5 2 — |C 2 In 3 2 + |£ 3 In 2 2 + 4 In 2Li 4 + 4Li, 5 Q) 

= -0.3851859504113720162670058305844... , 

LsLsc w (§) = = 0.6336780725876149549802789098028 ... , 

Ls L sc w (f) = -^vrC 2 Ls 2 (f) + HC 2 C 3 + f§§C5 + (f) 

-lb Ls 4 (f) - l Ls 2 (f) Ls 3 (t) + l Ls 5 1} (f) - 
= 0.7441484098381945153773320072924... , 

LsLsc 1j2)2 (f) = -SttC 2 Ls 2 (f) - |Ls 2 (|) Ls 3 (f) + J|7rLs 4 (f) 

~nl&A, (^) + - ^C5 - jigX5 

= -0.3438441037234077110928896695989... , 

LsLsc 1)2)2 (f) = ^7 tC 2 Ls 2 (f) + ±Ls 2 (f) Ls 3 (f) - (f) 

+ 24 7I 'Ls 4 (f) - HC 2 C 3 - WC5 - |Ls« (f) + £ Xb 
= -0.2676755241093223655003711724570... , 

LsLsco, 3 , 2 (j^) = — LsLsc 4i i )3 + ^7tC 2 Ls 2 + |Ls 2 Ls 3 + §||C 2 C 3 

—HC5 - ^7rLs 4 (f) - ivrLs 4 (f) + jLs^ (f) + 57^ 

+4Li5 + 4 In 2Li 4 — |^ 2 In 3 2 + |£ 3 In 2 2 + ^ In 5 2 

= -0.1728468405935728535297398287784... 
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LsLsc 0)3 ,2 (yy) = LsLsci,i )3 (§) - H7rC 2 Ls 2 (f) - |Ls 2 (f) Ls 3 + ^C 2 C3 
+ wCs + 04g7I"Ls 4 (f) - ^7 tLs 4 (f) - |Ls^ (f) - ±X5 
-12Li 5 (I) - 12 In 2Li 4 (|) + 2( 2 In 3 2 - f( 3 In 2 2 - § In 5 2 
= -0.2193285689662125608977695929575... 

LsLscj. 1 3 = — LsLsci^^ (f) + ^7rC 2 Ls 2 + |Ls 2 Ls 3 — IHC 2 C 3 

“lifCs - ^vrLs 4 (f) - ^vtLs 4 (f) + |Ls^ (f) + ± X 5 
+16Li 5 (i) + 16 In 2Li 4 (±) - |( 2 In 3 2 + 7( 3 In 2 2 + ^ In 5 2 
= 0.4544220738685365169842148301889.... (B.20) 

where the high-precision numerical evaluation of the generalized log-sine functions was per¬ 
formed with the lsjk-program [39] and 

LsLsci,i i3 (§) = 0.32403774485559909073259711644693958993 .... 
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